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Both the differential tunneling conductance and the surface local density of states (LDOS) of a 
[110]-oriented high-temperature superconductor in the pseudogap (PG) regime are studied theoret- 
ically. As a competing candidate for the mechanism of PG state, the charge-density wave (CDW), 
spin-density wave (SDW), d-density wave (DDW), and a!- wave superconducting (DSC) orderings 
show distinct features in the tunneling conductance. For the CDW, SDW, and DSC orderings, 
the tunneling conductance approaches the surface LDOS as the barrier potential is increased. For 
the DDW ordering, we show for the first time that there exist midgap states at the [110] surface, 
manifesting themselves as a sharp zero-energy peak in the LDOS, as in the case of DSC ordering. 
However, due to the particle-hole pair nature of the DDW state, these states do not carry current, 
and consequently the one-to-one correspondence between the tunneling conductance and the surface 
LDOS is absent. 
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One of the most intriguing properties of high tempera- 
ture superconductors is that, in the underdoped regime, 
the density of low-lying excitations is suppressed be- 
low a characteristic temperature T* , which could be 
well above the superconducting transition temperature 
T c |Q. The origin for this so called pseudogap (PG) 
phenomenon remains unsettled. The existing scenarios 
for the mechanism of PG state include: the preformed 
pair model, where the Cooper pairs are formed above 
T c but either the phase associated with them fluctuates 
so strongly Q or the coherence length (essentially the 
size of a Cooper pair) is so small || that the super- 
conducting transition temperature is reduced from T* 
to T c ; stripes @, 0, |[ or antiferromagnetic fluc- 
tuations []lo|| , relying on the spin and/or charge fluctua- 
tions; spin-charge separation flT|| , where the PG comes 
from the spinon sector; quantum criticality p2| , of which 
the d-density wave (DDW) ordering breaking both time- 
reversal and translational symmetry |l3[ and the circulat- 
ing currents breaking time-reversal but preserving trans- 
lational symmetry |l4j , are the most representative mod- 
els. Recently, several tunneling experiments have been 
proposed for both the preformed pa ir [Ji~5| , [l6| and anti- 
ferromagnetic correlation scenarios Jl7[. In this paper, 
based on a lattice model, we calculate the tunneling con- 
ductance, within the framework of a scattering approach, 
of a [110]-oriented normal- metal/undcrdopcd high-r c su- 
perconductor junction. We typically consider the pos- 
sible orderings: (i) the charge-density wave (CDW) to 
model the charge stripes; (ii) the spin-density wave 
(SDW) to model the spin stripes or antiferromagnetic 
correlations; (iii) the DDW state; (iv) the d-wave su- 
perconducting (DSC) state. The exposure of distinct 
features among these orderings may help us pin down 
the mechanism for the PG. On the other hand, it has 
been predicted that a sizable areal density of midgap 
states (i.e., states with essentially zero energy relative to 



the Fermi surface) exists on the [110]-oriented surface of 
a dj.2_ y 2-wave superconductor. These states gives rise to 
a narrow surface local density of states (LDOS) peak at 
the Fermi energy, where the bulk density of states dips 
to zero. One of the consequences of these midgap states 
is a sharp zero bias conductance peak in single particle 
tunneling [[l9| pc| ], which has been confirmed by several 
carefully controlled experiments |2^, [2^, Q . A ques- 
tion arises naturally: Whether such kind of midgap states 
also exist in those normal-state orderings, and if so, how 
to distinguish them from those in the DSC ordering. This 
question constitutes another motivation for the present 
work. In light of this, we also calculate the LDOS near 
the [110]-oriented surface of the superconductor with all 
orderings listed above. We choose the [110] orientation, 
where the crystalline a axis of the superconductor has 
a j angle with respect to the direction normal to the 
surface, because the surface sensitivity of the electronic 
structure of some (e.g., the DDW and DSC) orderings 
can be manifested most pronouncedly in this setup. 

The equation of motion for quasiparticles in a super- 
conductor with a variety of orderings can be described 
in a unified manner by the generalized Bogoliubov-de 
Gennes equations defined in a two-dimensional tight- 
binding model: 
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Here the quasiparticle wavefunction 



Ui,a 

Vi.-rr 



comprises 



the electrons at energy E with spin a coupled with 
the holes of the same energy but opposite spin —a; 
i = (ia,ib) are the lattice indices in the coordinate sys- 
tem defined by the crystalline a and b axes (Hereafter the 
lattice constant is chosen to be unity); the single parti- 
cle Hamiltonian Hij = —tSi +T .j + W t j + (Ui + Si,- a )Sij, 
where t is the nearest-neighbor hopping integral and 
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t = (T a ,Tb) = (±1,0) or (0, ±1) represent the relative 
position of those sites nearest neighboring to the ith site; 
W,j = (-l) ?; " +i6 VM 4i for the conventional CDW order- 
ing while Wij = {-if^iWdSgnirl - t 2 ) for the DDW 
ordering, Si )C = a{— l) la+lb So represents the SDW order 
parameter; U% is the scattering potential to model the im- 
purities, surface or interface. The quantity A^j = A s Sij 
for the s-wave superconducting order parameter while 

c 2 
i j 



Aij = A d sgn(x 2 



yfj) for the DSC order parameter. 



We note that the CDW, SDW, and DDW fields partic- 
ipate in the diagonal elements of the 2x2 BdG ma- 
trix equations in spin space while the superconducting 
order parameter appears as off-diagonal elements. Of 
most interest to us in the present work, as explained 
above, we consider the tunnel junction comprising a nor- 
mal metal on the left hand side and a superconductor 
on right hand side with a [110] oriented surface. We 
choose that the interface runs along the y direction while 
the normal to the interface is parallel to the x direction, 
along which the transport takes place. The insulating 
barrier at the interface is modeled by a scattering poten- 
tial, Ui = Uo5i x fi. In addition, we assume that the order 
parameter is identically zero in the normal metal and 
constant in the superconductor. In general, due to the 
proximity effect, the order parameter is degraded when 
the interface is approached from the interface within a 
characteristic length. However, the discontinuity of the 
pair field remains at the interface and the physics dis- 
cussed here is not changed qualitatively. By using the 
translation invariance of the system along the specular 
interface with spacing the differential conductance 
(in units of the conductance quantum e 2 /h) is calculated 
within the framework of the Blonder-Tinkham-Klapwijk 
(BTK) pj| theory: 
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where Nyy/2 is the linear dimension of the junction along 
the y direction, r n and r a are respectively the normal and 
Andrccv reflection amplitude for a beam of electrons with 
spin a incident from the normal metal at a fixed energy 
E and the transverse wavevector k y = y/2nl y /N y with 
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1. We determine the values of r n and 



r a by matching the wavefunction. Away from the in- 
terface, the solution to the BdG equations consists of the 
incoming electron, reflected and transmitted electron and 
hole waves. The normalized zero-temperature tunneling 
conductance is defined as: g(E) = G(E)/G nn (E), where 

G nn (E) = J^J2k y ,<T 4+4/«l With tA - = - 2tcOS (K/V2) 

is the tunneling conductance through a normal junction 
of the identical interface condition. We remark that 
the BTK formula can be considered as an extended ver- 
sion of the Landauer-Biittiker formula for normal 
metal/superconductor junctions. The latter has found 
spectacular success in the description of quantum cohcr- 
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FIG. 1: The normalized tunneling conductance (a) of a [110]- 
oriented N/CDW junction for various strength of the tunnel- 
ing barrier, Uo = lOt (red line), At (blue line), and It (green 
line). Also shown in (b) the LDOS near the [110]-oriented 
surface. The amplitude of the CDW, W = O.lt. 



ent transport in mesoscopic normal conductors. At the 
same time, we exactly diagonalize Eq. (|jj) to calculate the 
LDOS for a [110]-oriented surface of the superconductor, 
which is defined as: 
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where for a fixed value of k y , the summation is over the 
complete set of eigenstates. 

Tunneling into CDW. In the case of normal- 
metal/CDW junction, the reflection amplitudes are ob- 
tained as: 



(Uo 



i(E-Wo) > 



isgn(E) 



i(E-W ) 



isgn(E) 



(4a) 
(4b) 



Notice that for the CDW, SDW, and DDW orderings (as 
considered below), the electron correlation is through the 
particle-hole channel so that the charge is conserved and 
the Andreev reflection process is absent. In Fig. 0(a), 
we show the tunneling conductance g(E) as a function 
of energy for various strength of the tunneling barrier. 
The surface LDOS p(E) is displayed in Fig. 0(b). As 
shown in Fig. 0, both the tunnel conductance and the 
LDOS vanish when the energy is below the energy gap 
Wq. As the tunnel barrier becomes more opaque, the 
profile of the energy dependence of the tunneling con- 
ductance resembles that of the surface LDOS. Both g(E) 
and p(E) differ dramatically from those of a normal- 
metal/s-wave superconductor junction [ p5[ , by exhibit- 
ing no coherent peak at one of the gap edges. They are 
also asymmetric with respect to the Fermi energy E = 0. 
We also find the property that g(E, Wq) = g(—E, —Wo) 
and p(E, Wo) = p{—E, —Wo), which show the strong de- 
pendence of both the tunneling conductance and LDOS 
on the macroscopic phase of the CDW. 
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FIG. 2: The normalized tunneling conductance (a) of a [110]- 
oriented N/SDW junction for various strength of the tunnel- 
ing barrier, Uo = lOt (red line), At (blue line), and It (green 
line). Also shown in (b) the LDOS near the [110]-oriented 
surface. The amplitude of the SDW, So = O.lt. 
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FIG. 3: The normalized tunneling conductance (a) of a [110]- 
oriented N/DDW junction for various strength of the tunnel- 
ing barrier, Uo = lOt (red line), At (blue line), and It (green 
line). Also shown in (b) the LDOS near the [110]-oriented 
surface. The amplitude of the DDW, W d = O.lt. 



Tunneling into SDW. In the case of normal- 
metal/SDW junction, we arrive at the reflection ampli- 
tudes as: 
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In contrary to the CDW case, where the spin is de- 
generate for the tunnel conductance, the reflection am- 
plitude for the SDW is spin dependent. In partic- 
ular, there holds the following symmetry property of 
the tunneling conductance between different spin bands: 
g a (E,So) = g- a (— E 1 So). It follows immediately that 
the total tunneling conductance is symmetric with re- 
spect to the Fermi energy, i.e., g(E, So) = g(—E, So). As 
shown in Fig. ||, both the tunneling conductance and the 
LDOS vanishes when the energy below the energy gap So- 
The profile of g(E) also resembles to that of p(E) as the 
strength of the tunnel barrier is increased. In addition, 
the coherent peak shows up at the gap edge E = ±So 
in both g(E) and p(E). Therefore, the overall behavior 
of the tunneling conductance of the N/SDW junction is 
qualitatively similar to the normal-metal/s-wave super- 
conductor junction. 

Tunneling into DDW. In the case of normal- 
metal/DDW junction, the reflection amplitudes are ob- 
tained as: 
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where Wk = — 4W<2 sin(fc y /v / 2). The energy gap for both 
the CDW and SDW ordering, as discussed above, is con- 
stant on the Fermi surface in the momentum space. The 
bulk energy gap of the DDW ordering is momentum- 
dependent, and of the form cos k x — cos k y , which shows 



clearly that the DDW energy is closed along the nodal 
directions on the Fermi surface. In Fig. |^, we show the 
tunneling conductance of the N/DDW junction and the 
surface LDOS of the superconductor at the DDW state. 
Noticeably, a sharp peak appears at E = in the sur- 
face LDOS, which indicates the existence of the midgap 
states at the [110]-oriented surface of the DDW ordering. 
However, the tunnel conductance exhibits a V-shaped 
feature, which persists even when the tunnel barrier is 
opaque. The V-shape originates from the momentum 
dependence of the energy gap of the DDW ordering. In 
particular, the tunneling conductance does not show a 
zero-bias anomaly, it instead dips to zero. It can also 
be seen clearly that when E = 0, |r„| 2 = 1, which gives 
g(E) = 0, regardless of the strength of the tunnel bar- 
rier. Therefore, in the tunneling limit, the profile of 
the conductance is strikingly different from the surface 
LDOS. Physically, the energy gap of the DDW ordering 
comes from the particle-hole pairs. In contrast to the 
superconducting state, there exists no Andreev reflection 
process to convert, at the interface, the single particle 
current into the supercurrent through the condensate, so 
that all quasiparticles with energy below the momentum- 
dependent DDW gap are completely reflected. Conse- 
quently, although there exist midgap states at the [110]- 
oriented surface of the superconductor in the DDW or- 
dering, these states do not carry charge. However, as in 
the DSC case ]27|] , the midgap states in the DDW state 
can also lead to a giant magnetic moment when a strong 
in-plane magnetic field is applied, the study of which is 
beyond the scope of current work. At this stage, we can 
conclude that no midgap states exist at the surface of the 
superconductor with either CDW or SDW ordering, but 
they really show up in the DDW ordering. 

Tunneling into DSC. The tunneling conductance of 
normal-metal/DSC junctions, has been studied very in- 
tensively in the continuum model. Considering that 
the Andreev reflection process should take place even 
above T c in the preformed pair scenario, or in the spin- 
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FIG. 4: The normalized tunneling conductance (a) of a [110]- 
oriented N/DSC junction for various strength of the tunneling 
barrier, Uq = lOt (red line), 4t (blue line), and It (green line). 
Also shown in (b) the LDOS near the [110]-oriented surface. 
The amplitude of the DSC, A d = O.lt. 



charge separation picture (where the spinons form d-wave 
resonant- valence-bond state) , we believe it is still helpful 
to present, within the lattice model, the results for the 
DSC case. We arrive at the reflection amplitudes: 
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where A k = — 4Ad sm(k y /V2). As shown in Fig. ||, 
the conductance at zero energy does not vanish. In par- 
ticular, when the strength of the tunneling barrier is in- 
creased, a sharp peak occurs at zero energy in the conduc- 
tance. In the opaque limit of the tunneling barrier, the 
profile of the conductance resembles that of the surface 
LDOS, which also exhibits a zero-energy peak, indicating 
the existence of midgap states. As we have emphasized 
above, the superconducting energy gap originates from 
the formation of Cooper pairs due to the electronic cor- 
relation in the particle-particle channel. Therefore, the 
charge is not conserved in a superconductor, and the An- 
dreev reflection can happen in the interface between the 
normal metal and superconductor. It can be found from 
Eq. (0) that, at E = 0, r n , CT = and r a , a = 1, regard- 
less of the strength of the tunneling barrier. Therefore, 
the existence of the midgap states in the superconduct- 
ing state corresponds to the perfect Andreev reflection 
process. In the phase fluctuation scenario, the long-range 
phase coherence is absent but the short-range correlation 
is still retained, the localized nature of the midgap states 
should be robust against the phase fluctuations. There- 
fore, if the phase fluctuation is the mechanism for the 
PG, we expect that the zero-bias conductance anomaly 
is still observable. 

In summary, we have studied the tunneling conduc- 
tance and the surface LDOS of a [110]-oriented high-T c 
superconductor in the undcrdoped regime. The distinct 



features among the different kinds of ordering — CDW, 
SDW, DDW, and DSC, are investigated. We have shown, 
in the opaque limit of the tunneling barrier, that: (i) 
for the CDW ordering, asymmetric tunneling conduc- 
tance and surface LDOS with vanishing intensity within 
the gap; (ii) for the SDW ordering, symmetric tunneling 
conductance and surface LDOS with vanishing intensity 
within the gap; (iii) for the DDW ordering, V-shaped 
tunneling conductance but zero-energy-peaked surface 
LDOS; (iv) for the DSC ordering, zero-energy peaked 
tunneling conductance and surface LDOS. These differ- 
ent features can serve to distinguish the competing sce- 
narios for the mechanism of the PG. We have also found 
the existence of surface midgap states in the DDW order- 
ing, which is, to best of our knowledge, the first example 
occurring in a gapped normal system. Also in this case, 
the correspondence between the tunneling conductance 
and the surface LDOS is broken down. 
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